Electron Transport and Hot Phonons in Carbon Nanotubes 
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We demonstrate the key role of phonon occupation in limiting the high-field ballistic transport 
in metallic carbon nanotubes. In particular, we provide a simple analytic formula for the electron 
transport scattering length, that we validate by accurate first principles calculations on (6,6) and 
(11,11) nanotubes. The comparison of our results with the scattering lengths fitted from experi- 
mental FV curves indicates the presence of a non-equilibrium optical phonon heating induced by 
electron transport. We predict an effective temperature for optical phonons of thousands Kelvin. 

PACS numbers: 73.63.-b, 73.63.Fg, 63.20.Kr, 73.22.-f 
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Single wall carbon nanotubes (SWNT) have unique 
properties that make them strong candidates for future 
electronic devices. They can act as one-dimensional 
quantum wires with ballistic electron transport 0, 0, Q • 
Due to the strong C-C bond, they can carry the high- 
est current density of any material before they break. 
This makes them the best candidates as interconnects in 
integrated circuits and for high performance field effect 
transistors. It is thus essential to understand what ulti- 
mately limits the high current performance of nanotube 
devices. High field transport measurements have shown 
that electron-phonon coupling limits the ballistic behav- 
ior 0,0>0I- At high bias (>0.2V) electrons scatter with 
optical phonons, while at low-bias (<0.2V) they scatter 
with acoustic phonons. Ballistic transport is possible up 
to few hundreds nanometers in low bias, but the electron 
mean free path significantly drops at high bias 0] • 
By constructing devices smaller than the high bias scat- 
tering length a boost in performance is achieved 0] • 

The Boltzmann theory is used in Refs. 0, 0] to fit 
the scattering length (l op ) from the measured I-V curves. 
It is striking that the analysis of several SWNTs, with 
1-3 nm diameters, gives a similar Im, ~ 10—15 nm for 
optical phonons back-scattering 0, 0] . Alternatively, 
the scattering lengths could be derived knowing (i) the 
electron-phonon-coupling (EPC), and (ii) the phonon oc- 
cupation during transport. Thus, if the EPC is inde- 
pendently known, the phonon occupation can be derived 
from the experimentally fitted scattering length. A reli- 
able determination of the EPC is essential to apply this 
procedure. 

Several attempts to determine the EPC using tight- 
binding (TB) have been reported SUES 13 El Hi 
The predicted EPCs are not reliable since they strongly 
depend on the used TB parametrization. E.g., for the 
graphene unit cell and the E 2ff phonon, the maximum 
of the EPC square at T \D^ ax \ 2 , has been calculated to 
be 1L-42, 35/58, 164 (eV/A) 2 in Refs. ®, [30, 
[lfl Q] and Q respectively. These numbers give an 
order of magnitude spread in the calculated scattering 
length, being this proportional to |Z?p a:E | 2 . On the other 



hand, the optical phonons EPC were never directly mea- 
sured. But, we recently showed that graphite EPCs can 
be directly determined from the experimental phonon- 
dispersions or from the Raman D peak dispersion [lfij . 
The measured EPCs are in excellent agreement with ac- 
curate density functional theory (DFT) calculations ^3 • 
This confirms the reliability of DFT for determining the 
EPCs in graphitic materials, such as SWNTs. 

In this Letter, we demonstrate the key role of phonon 
occupation in limiting the high-field transport in nan- 
otubes. To do so, we use DFT to compute with high 
accuracy the optical-phonons EPC in (6,6) and (11,11) 
SWNTs. We demonstrate that curvature does not af- 
fect the EPCs in SWNTs with diameters >0-8 nm^uch 
as those used in transport measurements 0, 0, 0, 0] ■ 
Then, we prove that the EPCs for SWNTs of arbitrary 
chirality can be derived from those of graphite, by using a 
simple analytic formula. We use our computed and mea- 
sured 0| EPCs to obtain the electron mean free path 
for optical phonon scattering in high-field quasiballistic 
transport. We obtain the phonon occupation by compar- 
ing our scattering length with the / op fitted in transport 
measurements. 

SWNTs are identified by the chiral indexes (n,m), giv- 
ing the chiral vector Ch = nai + ma.2, with ai and &2 
the graphene lattice vectors 16]. A SWNT is obtained 
by folding graphene so that that the two atoms connected 
by Ch coincide. In real space, a SWNT is periodic along 
its axis according to the translational vector T [lfjj . For 
diameters >0.8 nm the electron states are very well de- 
scribed by the folded graphene bands ^3 ■ Thus, we can 
safely use zone-folding to describe the states relevant to 
charge transport. The SWNT electron-states correspond 
to the graphene states having periodicity Ch, i.e. such 
that k • C h = 2-7TO, o being an integer. The SWNT elec- 
tron states are labeled by two-dimensional momentum 
vectors k, which cut with a series of parallel lines the 
graphene Brillouin zone (BZ) [lfjj ]. 

The SWNT optical phonon frequencies are modified 
by the quantum confinement of electron states along the 
circumference For diameters >0.8 nm, at 



room temperature, the deviation of the SWNT phonon- 
frequencies from folded graphene is <5% |l9(. Such 
deviation has a negligible impact on scattering lengths. 
Thus, we use zone-folding to describe SWNT phonon fre- 
quencies and eigenmodes, and we label the phonons by a 
2-dimensional momentum vector q in the graphene BZ. 
The SWNT electron and phonon states are normalized 
in the 1-dimensional SWNT unit cell, with period T. 

The decay time for an electron state k, band i and 
energy e^i into the electron band j and energy e(k+qVp 
due to scattering with q or — q phonons in branch 77, and 
with energy fiw qj) , is given by the Fermi golden rule [20| : 
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where N is the number of q points, the two Dirac 5 distri- 
butions describe the processes of phonon absorption and 
emission, and n qr/ is the phonon occupation number. For 
thermal equilibrium, n q7) is the Bose- Einstein occupation 
factor riqr, = [exp(fiw q?) /fc B T) - Within DFT, the 

EPC is defined as #(k+q)j,ki = £>(k+q).7.ki-\/fi7 (2Mw qt) ), 
where M is the atomic mass, 
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|k, i) is the electron state, and AVq,, is the Kohn-Sham 
potential derivative with respect to phonon displacement. 

For optical phonons, the phonon dispersion is much 
smaller than the electron band dispersion. Thus, consid- 
ering only phonon emission, i.e. only the second S, Eq. ^ 
becomes 
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|£>(k+qb-,k t | 2 p[e(k+q)j] (n_ qJ7 + l), (3) 



where p[e(k+ q )j] is the electron density of states for band 
j, and q is fixed by e( k+q )j = e ki - fi^ qj) . 

For graphene, the electronic gap is zero for the n 
bands at the two equivalent BZ points K and K' = 2K. 
Metallic SWNTs have electron states corresponding to 
the graphene K and 2K points. For bias smaller than 
1.0 V, the mobile electrons belong to it bands near the 
Fermi energy e^- Both the initial and final states, k and 
k + q, are close to K or to 2K. The it bands can be 
considered linearly dispersive, crossing at ep with slopes 
13 and — /?, with (3 = 5.52 A-eV, and density of states 
p[ £ (k+q)j] = T/(2ir/3). Phonons connect electron states 
near ep only for q ~ Y of for q ~ K, thus, in Eq. J3J), 
ujq V can be approximated by lut v or by oj-k v - After the 
scattering, the electron can either maintain its propa- 
gation direction, forward scattering (fs), or reverse it, 
back scattering (bs). We thus have four possible scat- 
tering processes, Fig. ^ We label the corresponding de- 
cay times Tp s , Tjc, t|? , and t||- The total decay time 
is given by 1/r = l/r^ s + 1/t£ s + l/rf s + l/r|f. To 
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FIG. 1: Conduction band electrons in a metallic SWNT. First 
order decay induced by the emission of a phonon with fre- 
quency oj q . An electron state can be scattered into the same 
band by forward scattering (fs) or into another band by back 
scattering (bs). The resulting four possible processes are la- 
beled by the phonon vector (r or K) and by the type of 
scattering (fs or bs). ef is the Fermi energy. 



model the transport with Boltzmann 0, one needs to 
know the decay times for each individual process. Pos- 
sible scattering into higher energy bands is neglected. 
This assumption is valid for electron energies lower than 
2/3/d, d being the tube diameter. Note that the electron- 
phonon scattering processes involved in electron trans- 
port are entirely analogous to those involved in Raman 
spectroscopy [la I22I |23| . 

We now compute the SWNT EPCs by folding the cor- 
responding graphene EPCs. We then demonstrate the 
validity of this approach, which, so far, has never been 
proven. The EPC (D) for a SWNT of arbitrary chirality 
can be obtained from the graphene EPC (D), using: 
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where S = Oq-\/3/2 is the graphene unit-cell surface, ao = 
2.46 A is the graphene lattice spacing, and 5=|Ch X T| 
is the nanotube unit-cell surface 0. Eq.H comes from 
Eq.[21 using the appropriate normalization of phonon and 
electron states. Combining all the previous results, for 
metallic SWNTs, Eq.|3 becomes: 
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where d — C^/ir — aov n 2 + nm + m 2 /it is the diameter. 

For electron states near ep in graphene, |Z)P_is not 
negligible only for the Y — E 2s and K — A' : modes [lj|- For 
optical phonons near T, the doubly degenerate Yi2 g mode 
splits into two almost longitudinal (LO) and transverse 
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TABLE I: \D\ 2 between tt states around K. 
(TO) modes. For small k', the corresponding EPCs are: 

= (5£) P [i±cos(9 + e r )} 

= (D 2 ) F [l T cos(0 + 0')](6) 



-,LO/TO |2 
(K+k'+q)7r*,(K+k')7rl 

|fjLO/TO ,2 
l- t ^(K+k'+q)7r*,(K+k')7r* 1 



with 9 the angle between k' and q, 9' the angle between 
k' + q and q, and tt and tt* are the lower and higher tt 
bands ^J. + or — refer to the LO and the TO modes, 
respectively. For the K — mode we have: 



|-D(2K+k'+q')7r*,(K+k')ir| 2 



\D 



(2K+k'+q')*-*,(K+k')ir* 



(Dfov (l + cosfl") 
(^k)f (1-cosO.(7) 



where 9" is the angle between k' and k' + q' |24( . From 
our DFT calculations, = 45.60 (eV/A) 2 , and 

(D^) F = 92.05 (eV/A) 2 Q. These values reproduce 
very accurately the EPCs extracted from the measured 
graphite phonon dispersion 15]. 

From Eas. l6l7l we obtain the EPCs for all the scattering 
processes of Fig. ^ The results are in Tab. D] Note that 
for metallic SWNTs, k' and q or q' are always parallel 
to the axis, and 9, 9' and 9" are always equal to or tt. 
Thus, the results of Tab. [U do not depend on chirality. 

Using Eq.[S]and Tab.[I]we can easily compute the scat- 
tering length I (I — t/3/H) for a generic metallic SWNT. 
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In Refs [3- 0. 0] phonons are assumed thermalized at 
room temperature (n-q„ ~ 0) and the contribution of 
forward scattering is neglected. Thus, to compare with 
experiments, we consider only back scattering and define 
lop = (1/^k + Vr-Lo) 1 - Using the graphene phonon 
frequencies (wr = 196.0 meV, wk = 161.2 meV) [l5j . 
we obtain the results of Tab. [H] In particular, we get a 
simple scaling between scattering length and diameter: 



l op = 65 d. 



(9) 



To validate these results, we perform an explicit cal- 
culation of \D\ 2 on (6,6) and (11,11) armchair SWNTs. 
DFT calculations are done with the gradient corrected 
functional of Ref. . We use plane- waves (40 Ry cut- 
off) and pseudopotential approaches [i^] . Electron states 
are occupied with an Hermite-Gauss smearing of order 1 
and a width of 0.05 Ry [2?J. EPC calculations are done 



TABLE II: Scattering lengths (in nm) for the processes in first 
column. Phonon occupation is assumed thermalized at room 
temperature. All lengths are in nm. In parenthesis, results 
from direct DFT calculations on (6,6) and (11,11) SWNTs 
with diameter d = 0.8 and 1.49 nm. The other values are 
obtained by folding the graphene EPCs. 



using the perturbative method of Ref. |2j| . We use an 
hexagonal super-cell with a neighboring tube distance of 
5 A and a grid of 16 k-points with a logarithmic distribu- 
tion around the K point. In Tab. ITU we report the scatter- 
ing lengths obtained using the \D\ 2 computed explicitly 
for the two tubes. Tab. [Tj] shows that the zone-folding 
results are quite accurate for the 1.5 nm diameter tube, 
and only a <20% difference is present for the smaller- 
diameter tube, due to curvature effects. 

Thus, Eq. IS] combined with our computed and mea- 
sured graphite EPCs [l^, provides an independent as- 
sessment of the scattering lengths for the 1-3 nm diam- 
eter SWNTs of Refs. However, the scattering 
length reported in Refs. 0,|a,E| is ~10-15 nm. This is 
one order of magnitude smaller than what is predicted by 
Eq.|Hl This disagreement cannot entirely be attributed to 
the experimental uncertainty in measuring the scattering 
lengths, since three independent experiments on a variety 

^^Ti 

of different tubes obtained very similar results |4j, |5|, |6j . 
Also, it cannot be ascribed to an error in the DFT EPCs, 
since our computed EPCs reproduce very well those ex- 
tracted from the experimental graphite phonon disper- 
sions or from the Raman D-peak dispersion Thus, 
we conclude that the hypothesis of thermalized phonon 
occupation (n_ q?) ~ in Eq. |HJ does not hold. 

A significant phonon occupation n can explain the 
small value of the measured scattering length l op . This 
hot phonon generation can occur if, during high-field 
electron transport, the optical phonons excitation rate 
is faster than their thcrmalization rate. With a high 
phonon occupation, both phonon emission and absorp- 
tion processes are equally relevant. We can estimate the 
scattering lengths by assuming that the phonon occu- 
pation n is independent of q and 77, and by summing 
the absorption and emission contributions. The scatter- 
ing lengths are then obtained by substituting (n_ qr; + 1) 
with (2n + l) in Eq. El 



/ op = 65 d/(2n + l). 



(10) 



n in the 2.7-5 range is necessary to reconcile the scat- 
tering lengths derived from the computed and measured 
EPCs and form the fit of the measured IV curves. This 
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corresponds to an effective temperature for the occu- 
pation of optical modes > 6000 K. This temperature 
is related only to the phonons directly excited by first 
order scattering, and not to other phonons, otherwise 
the SWNT would melt. Thus, during high-bias electron 
transport, the phonons are not in thermal equilibrium. 

These results are consistent with the observation 
that high-bias saturation currents in SWNTs on a sub- 
strate are significantly higher than those in suspended 
SWNTs H^]. Indeed, the effective temperature of op- 
tical phonons in suspended SWNTs is expected to be 
higher due to the absence of a thermally conductive sub- 
strate for heat sinking. Also, it was recently reported 
that electrons inelastically tunnelling from a scanning 
tunnelling electron microscope tip into a SWNT induce a 
non-equilibrium occupation of the radial breathing mode 
phonons |30j . In this case, the tunnelling electrons are 
mostly coupled with the radial breathing mode, since the 
associated vibrations directly modulate the tip-nanotube 
distance. On the contrary, in high-field transport exper- 
iments, electrons mostly excite optical phonons, which 
are the dominant scattering source. 

The generation of hot optical phonons, during high 
field electron transport, can be directly detected by a Ra- 
man scattering experiment combined with electron trans- 
port. The occupation of a given phonon mode can be de- 
termined from the ratio of the corresponding Stokes and 
anti-Stokes peak intensities. A rise in this ratio would 
be direct evidence and would allow to quantify the hot 
phonons generation. 

Finally, our approach can be easily extended to cal- 
culate scattering times and mobility in semiconducting 
SWNTs. In this case, the scattering times depend on 
electron energy and chirality since the bands are hyper- 
bolic and the electron density of states is energy depen- 
dent. The angles 9, 9', and 9", entering in Eqs. ® and 
0, are not equal to or ir, as in metallic SWNTs, but 
depend on the chirality and on the initial and final state 
momenta (k and k + q) . These angles can be computed 
by simple geometrical considerations. 

In conclusion, accurate DFT calculations of EPCs in 
graphene, (6,6) and (11,11) SWNTs, combined with a 
simple zone-folding model and the experimental graphite 
EPCs are used to interpret the saturation SWNTs I-V 
curves of Refs. We show that the optical phonon 

occupation is greatly increased during high-bias electron 
transport, with an effective temperature of thousands K. 
Such a high temperature causes a strong reduction of the 
ballistic scattering length. This suggests coupling the op- 
tical phonon mode with an heat sink in order to reduce 
their effective temperature. This would increase the scat- 
tering length up to the maximum value of EqO which 
sets the ultimate limit of ballistic transport. Finally, the 
phonon generation by electron scattering is analogous to 
the photon generation by stimulated emission in semicon- 
ducting lasers. This suggests that SWNTs, under high 



bias, can act as a possible source of coherent phonons. 
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